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Abstract 

We show that a combined analysis of B — ► Kvv and B — > Kl + l~ allows for new 
physics tests practically free of form factor uncertainties. Residual theory errors 
are at the level of several percent. Our study underlines the excellent motivation 
for measuring these modes at a Super Flavour Factory. 
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1 Introduction 



The standard model of quark flavour physics has successfully passed all experimental 
tests to date. This includes the observation of a substantial number of rare processes and 
CP asymmetries, which are consistently accounted for within the Cabibbo-Kobayashi- 
Maskawa (CKM) description of quark mixing. On the other hand, many essential fea- 
tures of the standard model, most notably in the flavour sector, are still not satisfactorily 
understood on a more fundamental level. Deviations from standard expectations, which 
could guide us towards a better understanding, appear to be small in general in view of 
the basic agreement between theory and observations. In this situation precision tests of 
flavour physics become increasingly important, which motivates current efforts to build a 
Super Flavour Factory [1,2]. Such a facility will enable an exciting program in B physics 
[3,4]. 

One of the best opportunities in this respect could be provided by the study of 
b — > svv transitions, induced by interactions at very short distances. Theoretically 
ideal would be an inclusive measurement of B — > X s vv, where the hadronic matrix 
element can be accurately computed using the heavy-quark expansion. Unfortunately, 
because of the missing neutrinos, an inclusive experimental determination of the decay 
rate is probably unfeasible. More promising is the measurement of exclusive channels 
such as B — » Kvv, B — > K*vv. In this case a clean theoretical interpretation requires, 
however, the control of nonperturbative hadronic form factors. Direct calculations of 
form factors suffer from sizable uncertainties. Additional experimental input to eliminate 
nonperturbative quantities can therefore be very useful. Important examples are the 
related decays K — > ttvv, where the hadronic matrix element can be eliminated with the 
help of K + — > 7r°e + z/ using isospin symmetry. As discussed in [5], a similar role could be 
played by the semileptonic mode B — > irev for the rare decay B — > Kvv. This strategy 
is limited by the breaking of SU (3) flavour symmetry of the strong interaction, which is 
also difficult to estimate with high accuracy. 

In this paper we propose to perform a combined analysis of the rare decays B Kvv 
and B — > Kl + U . As we shall discuss, this option has several advantages for controlling 
hadronic uncertainties. It allows us to construct precision observables for testing the 
standard model and for investigating new physics effects. In particular neither isospin 
nor SU (3) flavour symmetry are required and form factor uncertainties can be eliminated 
to a large extent. 

The paper is organized as follows. Section 2 summarizes the experimental status. Sec- 
tion 3 collects basic theoretical results. It includes a discussion of B — > K form factors, 
weak annihilation and nonperturbative corrections in B — > Kl + U , and the background 
for B~ — > K~vv from B~ — > t~v t — > K~v T v T . Precision observables are discussed in 
section 4. Section 5 comments on the effects of new physics and conclusions are pre- 
sented in section 6. Further details on form factors relations and on weak annihilation 
are collected in the appendix. 
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Table 1: The normalized g 2 -spectrum for B — > Kl + l~ . Shown are the partial 
branching fractions in six bins of q 2 (or s = q 2 /m 2 B ) from [10], normalized by 
the central value of the integrated branching fraction in (3). These quantities 
are denoted by AB/B in the table. 



g 2 [GeV 2 ] 


s 


AB/B 


0.00-2.00 


0.00-0.07 


0.169 ±0.038 


2.00-4.30 


0.07-0.15 


0.096 ±0.027 


4.30-8.68 


0.15-0.31 


0.208 ± 0.042 


10.09-12.86 


0.36-0.46 


0.115 ±0.031 


14.18-16.00 


0.51-0.57 


0.079 ±0.033 


> 16.00 


> 0.57 


0.204 ±0.042 



2 Experimental status 

In this section we summarize briefly the current experimental situation. For the branch- 
ing ratios of the neutrino modes B — > Kuu only upper limits are available at present. 
They read [6,7,8,9] 

B(B- -> K~vv) < 14 • 10~ 6 (1) 
B(B° -> K°uu) < 160 • 10~ 6 (2) 

Here CP averaged branching fractions are understood. We note that the limit is more 
stringent for the B~ channel. 

The most accurate experimental results for B — > Kl + l~ are from Belle [10]. The 
extrapolated, non-resonant branching fraction is measured to be 

B(B -> Kl + T) = (0.48^ ± 0.03) • 10~ 6 (3) 

consistent with results from BaBar [11]. The recent paper [10] also contains information 
on the g 2 -spectrum in terms of partial branching fractions for six separate bins. The 
results for the normalized g 2 -spectrum, adapted from [10], are given in Table 1. 

3 Theory of B — > Kvv and B — ► Kl + l~ 

3.1 Dilepton-mass spectra and short-distance coefficients 

We define the kinematic quantities 

2 2 
1 m K , A \ 

s = — r K = —^- 4 



where q 2 is the dilepton invariant mass squared and is the mass of the B meson. 
The kinematical range of q 2 and its relation with the kaon energy Ex are given by 



Am 2 ~ < q 2 < {m B — m K ) 2 q 2 = m 2 B + m 2 K — 2rriBE K (5) 

We also use the phase-space function 

X K (s) = l + r 2 K + s 2 -2r K -2s-2r K s (6) 

The differential branching fractions for B — > Kvv and B — > Kl + l~ can then be 
written as follows: 

dB ^2 R ^ = r B ?^\VM^T( S )fl( S ) \<K m) ? (7) 

dB[B ^ Kl+n = r B ^^\V t ,V^\f (s)f + (s) (\a,(KU)f + Ko(A'«)H (8) 

Here tb is the £>-meson lifetime, Gf the Fermi constant, a = 1/129 the electromagnetic 
coupling and Vt s , Vtb are elements of the CKM matrix. A second contribution to the 
amplitudes proportional to V* s V u b has been neglected. It is below 2% for B — > Kl + l~ 
and much smaller still for B — > Kvv. 

The factorization coefficient a(Kvv) is simply given by a short- distance Wilson co- 
efficient at the weak scale, C L , [5] 

a(Kvv) = C V L = -—L- Vx X (x t ) (9) 
sin 6 vi/ 

where Xo is an Inami-Lim function [12] and x t = mf/M^, with m t = fht{mt) the 
MS mass of the top quark. The factor r]x = 0.994 accounts for the effect of 0(a s ) 
corrections [13]. At this order the residual QCD uncertainty is at the level of 1-2% and 
thus practically negligible. 

The factorization coefficient ag(KU) contains the Wilson coefficient Cg(fi) combined 
with the short- distance kernels of the B — > Kl + l~ matrix elements of four-quark op- 
erators evaluated at fi = 0{mi > ). The coefficient ag(Kll) multiplies the local operator 
(sb)v-A(Jl)v- At next-to-leading order (NLO) the result can be extracted from the ex- 
pressions for the inclusive decay B — > X s l + l~ given in [12,14,15], where also the Wilson 
coefficients and operators of the effective Hamiltonian and further details can be found. 
The NLO coefficient reads 

a 9 (Kll) =C 9 + h(z, s) (Ci + 3C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) 

~h(l,s) (4C3 + 4C4 + 3C5 + G0 

~h(0, s) (C 3 + 3C 4 ) + I (3C 3 + C 4 + 3C 5 + C 6 ) + ^C 7 (10) 
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Here 

CM = P + - 4Zo(x t ) + P E E (x t ) (11) 

sin k3 w 

is the Wilson coefficient in the NDR scheme, Pq, Pe are QCD factors and Eq, Yq, Zq are 
Inami-Lim functions. The function h(z, s), z — m c /mb, s = q 2 /m 2 arises from one-loop 
electromagnetic penguin diagrams, which determine the matrix elements of four-quark 
operators. In contrast to Cg the quantity a 9 (KU) is scale and scheme independent 
at NLO. To this order the coefficients Cj, % — 1, ... 7 in (10) are needed only in leading 
logarithmic approximation (LO). Note that here the labeling of C\ and C<i is interchanged 
with respect to the convention of [12]. 
The coefficient a w (Kll) is 

a 10 (Kll) = C w = - , I Y (xt) (12) 
sin cV 

3.2 Form factors 

The long-distance hadronic dynamics of B — > Kvv and B — > Kl + l~ is contained in the 
matrix elements 

(K(p')\sYb\B(p)) = U(s) (p + p'r + [fo(s) - f + (s)] ml ~f K q» (13) 

q 2 

(K(p')\sa^b\B(p)) = i fT [ s) [(p + pTq u -q^P + p'T] (14) 

which are parametrized by the form factors /+, fo and fx- Here q = p—p' and s = q 2 /m B . 
The term proportional to g M in (13), and hence fo, drops out when the small lepton 
masses are neglected as has been done in (7) and (8). The ratio fr/f+ is independent 
of unknown hadronic quantities in the small-s region due to the relations between form 
factors that hold in the limit of large kaon energy [16,17] 

Ms) m B + m K 

-r-TT = h 0{a s , A m b ) (15) 

/+(«) m B 

Here we have kept the kinematical dependence on tuk in the asymptotic result. In 
contrast to / + the form factor fx is scale and scheme dependent. This dependence is 
of order a s and has been neglected in (15). Within the approximation we are using we 
may take \x = rrib to be the nominal scale of fr. 

We remark that the same result for fx / /+ is also obtained in the opposite limit where 
the final state kaon is soft, that is in the region of large s = 0(1). This follows from 
the asymptotic expressions for /+ and fx in heavy hadron chiral perturbation theory 
[18,19,20,21,22]. From this observation we expect (15) to be a reasonable approximation 
in the entire physical domain. This is indeed borne out by a detailed analysis of QCD 
sum rules on the light cone [23], which cover a range in s from to 0.5. Relation (15) is 
further discussed in appendix A. 
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The ratio fr/f+ enters (10) as a prefactor of C7 from the matrix element of the 
corresponding magnetic-moment type operator Q 7 [5,12]. In writing (10) the relation 
(15) has already been used to eliminate fr/ f+- Since the C7 term contributes only about 
13% to the amplitude ag(Kll), the impact of corrections to (15) will be greatly reduced. 
A 15% uncertainty, which may be expected for the approximate result (15), will only 
imply an uncertainty of 2% for ag(Kll) or the B — > Kl + l~ differential rate. In practice, 
this leaves us with the form factor f+(s) as the essential hadronic quantity for both 
B -> Kvv and B -»• Kl + l~. 

The main emphasis of the present study is on the construction of clean observables, 
which are, as far as possible, independent of hadronic input. We will therefore consider 
suitable ratios of branching fractions where the form factor f+(s) is eliminated to a 
large extent. In order to assess the residual form factor uncertainties in these cases, 
but also to estimate absolute branching fractions, it will be useful to have an explicit 
parametrization of the form factor at hand. We employ the parametrization proposed 
by Becirevic and Kaidalov [27] in the form 

/+(s) s /+(0) \l Xki - M* = /+<0)|1 + tto6 ° S + ° {s2)] (16) 

The parameter b is given by 

2 

b = « 0.95 for m B * = 5.41 GeV (17) 

771 d* 

s 

b represents the position of the B* pole and will be treated as fixed, following [27]. The 
remaining three parameters do, b\ and /+(0) have been determined from QCD sum rules 
on the light cone (LCSR) in [23] 

/ + (0) = 0.304 ±0.042, OflWl.5, h = b (18) 

We will treat all three as variable parameters. This also includes 61, slightly generalizing 
the expressions from [23] where b\ is fixed at bo- The value for /+(0) in (18) is obtained 
from the relation [23] 

/ + (0) = 0.331 ± 0.041 + 0.25(ai(l GeV) - 0.17) (19) 

using the updated value [24,25] for the Gegenbauer coefficient a±(l GeV) = 0.06 ± 0.03 
as quoted in [26]. 

With bo fixed, the parameter do introduced in (16) determines the slope of the form 
factor at small s. We remark that the LCSR method is appropriate for the low-s region, 
which will be of particular interest for us. For completeness we give the relation of our 
parameters /+(0), do, b\ to the original parameters c B , ol = as, 7 = 7b from [27]: 

/+(0) = cb(1 - a B ) , a = ^- — , r = 7s (20) 
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As discussed in [27], the large-energy limit for the kaon implies the relation 7# = 1/olb 
or, equivalently, a b — b + b\. 

We determine next our default ranges for the shape parameters a and b\, which 
will be employed in the subsequent phenomenological analysis. Three main pieces of 
information will be used: The experimental data on the q 2 spectrum in Table 1, the 
LCSR results in (18), and asymptotic results for the form factor at maximum s, 



1-™X (21) 
m B J 



The third constraint will lead to a relation between a and b±. It follows from the 
asymptotic expression for f + (s m ) 

UM = 2^Ta) a = m ^ " mB (22) 

which can be derived within heavy-hadron chiral perturbation theory [18,19,20,21,22]. 
The largest uncertainty in (22) is due to the BB*K coupling g, sometimes also normalized 
as 9bb*k = 2ms5'//x, which is not known precisely. For the analogous, SU(3) related 
quantity gsB*^ a range of gBB*n — 42 ± 16 is quoted in [27]. This corresponds roughly 
to g — 0.6 ± 0.2. In view of this uncertainty, and since the main purpose here is the 
estimate of typical numbers, we have neglected subleading corrections to (22), which 
may be sizeable [20]. We recall that g is of order unity in the large-m_e limit. 
Equating (22) with f+(s m ) from (16) and using 

!_ W= a(!!!«L±A> (23) 



m B 



we obtain 



l-,„,(fr + ^-aA) = ^fe 2S±U) (24) 

1 - s m bi J+{0)JK 

We note that the denominator of the first term in (24), 1 — s m bi, scales as l/m^ in the 
heavy-quark limit, whereas the numerator remains of order unity (if is not too far 
below its typical value of 1.5). The first term then scales asm B , consistent with the 
heavy-quark scaling of the second expression. 

The constraint (24) can also be put in the form 

c «~ l = ^K < 25 > 

s m bo bo 

Numerically we have l/(s m b ) = 1.279. Within the uncertainty of c , displayed in (24) 
above, (25) implies a correlation between the shape parameters ao and bi/b . 

Independently of such theory constraints we might ultimately want to extract the 
form factor shape from experimental data. In this spirit, we have investigated how well 
different values of (ao,&i/&o) fit the current Belle measurements of the dilepton-mass 
spectrum in B — > Kl + l~ . For this purpose we show in Table 2 the x 2 -function 
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Table 2: Values of x 2 f° r various combinations of the form factor shape pa- 
rameters do and bi, determined from a comparison with the Belle data on the 
q 2 spectrum of B — > Kl + l~ (see Table 1). 



^"^^ a 
61/60 


1.0 


1.2 


1.4 


1.6 


1.8 


2.0 


0.5 


20.2 


14.7 


11.0 


8.8 


7.5 


7.1 


0.6 


20.2 


14.2 


10.3 


8.0 


6.9 


6.7 


0.7 


20.2 


13.5 


9.4 


7.1 


6.2 


6.4 


0.8 


20.2 


12.7 


8.3 


6.2 


5.7 


6.4 


0.9 


20.2 


11.8 


7.1 


5.3 


5.5 


7.0 


1.0 


20.2 


10.5 


5.7 


4.7 


6.2 


9.2 



X >(a M = ± {Vi - F ' ( :° M)2 (26) 

1=1 a% 

where the i — 1, . . . , 6, are the experimental values for the normalized, partial branch- 
ing fractions AB/B in each of the six bins, the <7j are the corresponding errors, and the 
Fi are the theoretical expressions depending on (cto, 61). It is clear that the experimental 
data are at present not accurate enough to allow for a precise determination of the form 
factor shape. However, the situation should improve in the future. At the moment our 
analysis merely serves to illustrate the general method. Nevertheless, some regions of pa- 
rameter space are already disfavoured, in particular low values of a , which, as discussed 
above, is consistent with theoretical expectations. We also observe that the best fit is 
obtained for a = 1.6, 61/60 = 1. A comparison of the best-fit shape of the theoretical 
spectrum with the Belle data is shown in Fig. 1 

Combining the information above, we adopt the following default ranges for the shape 
parameters 

1.4 < a < 1.8 0.5 < 61/60 < 1.0 (27) 

with 

a = 1.6 6i/6 = 1.0 (28) 

as our reference values. Within the range (27) for and 61/60 the parameter cq in (25) 
takes values between 1.5 and 3.9, compatible with (24). Our default parameters (28) are 
also consistent with the LCSR results [23] quoted in (18). 

3.3 B — ► Kl + l~: Weak annihilation 

As pointed out in [28], the exclusive decay B — > Kl + l~ receives contributions from weak 
annihilation diagrams already at leading order in the heavy-quark limit. In spite of this 
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(dB/ds)/B 




Figure 1: The shape of the B — > Kl + l~ spectrum, (dB/ds)/B, from the Belle data 
summarized in Table 1 (crosses) and from theory with the best-fit shape parameters 
a = 1.6, bi/b = 1 (solid curve). 



their impact is numerically small because of a strong CKM suppression (for the charged 
mode) or small Wilson coefficients (for the neutral mode). In this section we will quantify 
the size of weak annihilation. Since this contribution is small we work to leading order 
in a s . 

First we consider the case of B° — > K°l + l~ , and the region of low s ~ A/mj,. Weak 
annihilation can then only come from penguin operators giving rise to the transition 
bd — > sd. A virtual photon emitted from one of the quarks produces the lepton pair. 
The leading annihilation contribution in the heavy-quark limit is generated by the gluon- 
penguin operators Q3 and Q4 [12]. This effect can be evaluated using the methods of 
QCD factorization [28,29] as discussed in more detail in appendix B. The resulting 
correction to the coefficient ag(KU) reads 

Aa 9tWA ,34 = \C A + -C3 - / du) (29) 

V 3 J m B J + {s) J u - sm B -ie 

Here fs, Jk are the meson decay constants and Qd = —1/3 is the charge of the down- 
quark in the initial state. The leading light-cone distribution amplitudes of the B meson 
can be expressed by two functions, <P±{uj), of which only <p_ enters the integral in (29). 
For annihilation contributions scale dependent quantities such as the Wilson coefficients 
63,4 will be evaluated at a hard-collinear scale of fih = V where \i = 0(rrib) and 
A h = 0.5 GeV, following [29]. 

Our result (29) for weak annihilation agrees with eq. (68) of [28] , once it is adapted 
to the case of a pseudoscalar K meson as described immediately after eq. (69) of [28]. 1 



1 A factor of (— 2mb/Ms) is missing in front of the annihilation term on the r.h.s. of eq. (41) in 
We thank Thorsten Feldmann for confirmation. 
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For estimates of the annihilation effect we employ the model functions [30] 

= e -b>/»* = — e -"> M> (30) 

where o;o = O(Aqcd) serves to parametrize the uncertainty related to <p±. A summary 
of general properties of these wave functions can be found in [28,17]. They are satisfied 
by the parametrizations in (30). With (f>-(u) from (30), and denoting by Ei(z) the 
exponential integral function, the integral in (29) can be expressed as 



Ad 1 . (s) = I" duo = — e ~ smB '^ \-EUsm B /uj ) + in] (31) 

J u - s m B - it uo 

At a typical value of s = 0.1 we have ag = 3.96 + 0.05? and Aag t wA,34 = —0.036 + 0.034?. 
The correction (29) is seen to reduce the real part of ag by a small amount, which leads 
to a corresponding reduction of the branching fraction. This holds if s > 0.37uo/mB ~ 
0.025. Since the imaginary part of ag is much smaller than the real part, its impact 
on the decay rate is entirely negligible. Practically it is thus of no consequence that 
the imaginary part of the correction is comparable to the one of ag and that Imag is 
rather uncertain. Of particular interest for us is the size of the annihilation effect on the 
partially integrated branching fraction. For ujq = (0.350 ± 0.150) GeV, m^/2 < /i < 2rrib 
the reduction of the branching fraction integrated within 0.03 < s < 0.25 is at most 1%, 
which is indeed very small. 

Weak annihilation contributions to B° — > K°l + l~ also come from the remaining two 
QCD penguin operators Q5 and Qq. Because these have a chiral structure different from 
Q^i Qi) their contribution to weak annihilation is formally suppressed in A/m&. It turns 
out, however, that the suppression is not very effective numerically in this particular 
case. A similar situation is familiar from the factorizing matrix elements of Q 5 and Q 6 
for B decays into a pair of light pseudoscalar mesons [29]. The explicit calculation of 
the annihilation correction to a 9 (KU) from Q 5 and Q 6 proceeds as before and yields 



Aa w = - ( C 6 + ic.) lfa ^f/r* f ** *- M ■ (32) 

3 / m%f+{s) J uj-sniB-ie 

In comparison to (29) the correction in (32) carries, apart from the Wilson coefficients, 
a relative factor of — 2/i^/mB ~ —0.75, where \ik — — iTil/(m u + ma) [29]. Here 
m u , nid are the MS masses evaluated at the scale As anticipated, this relative 
factor is not small, although it is of 0{h./m^). A typical value for (32), at s = 0.1, is 
Aag^wA,56 — 0.045 — 0.043i. The sign of the real part is opposite to the case of Aag^A.34 
such that there is a tendency of the two contributions to cancel. Because |C 6 + C 5 /3\ 
is larger than |C 4 + C 3 /3|, it is possible that (32) even dominates over (29). Of course, 
(32) is formally a power correction and other power corrections to weak annihilation 
do exist. However, the numerically large factor 2^ik/^b is special to the annihilation 
matrix element of Qq = —2(db)s-p(sd)s+p with the presence of (pseudo) scalar currents. 
Taking into account the single power correction (32) thus appears justified. Adding 



9 



both corrections, (29) and (32), the net effect on the branching ratio integrated from 
s = 0.03 to 0.25 is a tiny enhancement. This enhancement remains below about 0.3% 
for lu = (0.350 ± 0.150) GeV and m b /2 < /i < 2m b . 

To summarize, the weak annihilation contributions to B° — > K°l + l~, which arise 
from QCD penguin operators, are negligibly small in practice, even though they are a 
leading-power effect. If the presumably dominant (chirally enhanced) power correction 
(32) is also included, the total impact of weak annihilation is further reduced. 

In the case of the charged mode B~ — > K~l + l~ the weak annihilation terms from 
QCD penguin operators are given by the expressions (29) and (32) with the replacement 
of the quark charge Qd — > Q u . Numerically the effect then receives an additional factor 
of —2, which still yields a negligible correction at the level of about one percent. Weak 
annihilation through the tree operators Q\ and Q 2 , which exists only for the charged 
mode, comes with large Wilson coefficients, but also with a strong Cabibbo suppression. 
This correction reads 

= ( Cl + u) !^!%^ r ^ — <l^l_ ,33) 

v ts v tb V 3 J m B f+{s) J u-sm B -it 

With (33) the branching ratio integrated from s = 0.03 to 0.25 is reduced by less than 
about 0.6% for u = (0.350 ± 0.150) GeV and m b /2 < fi < 2m b . Taking the three 
contributions from Qi,2, Qz,a and together, the reduction remains below 1%. Within 
an uncertainty of this order, weak annihilation is therefore negligible for B~ — > K~l + l~ 
as well. 



3.4 B — > Kl + l : Nonperturbative corrections 

In this section we comment on the theoretical framework for B — > Kl + l~ and on non- 
perturbative effects beyond those that are contained in the form factors. 

It is well known that, because of huge backgrounds from B — > Kip^ — > Kl + l~, the 
region of q 2 containing the two narrow charmonium states ip = ip(lS) and ip' = ip{2S) 
has to be removed by experimental cuts from the q 2 spectrum of B — > Kl + l~ . The 
overwhelming background from ip and vp' is related to a drastic failure of quark-hadron 
duality in the narrow-resonance region for the square of the charm-loop amplitude, as 
has been discussed in [31]. Nevertheless, the parts of the q 2 spectrum below and above 
the narrow-resonance region remain under theoretical control and are sensitive to the 
flavour physics at short distances. A key observation here is that the amplitude is largely 
dominated by the semileptonic operators 

Q 9 = (sb) V - A (ll)v 

Q w = {sb) v ^ A {ll) A (34) 

which have large coefficients Cg and Cio- These contributions are perturbatively calcula- 
ble up to the long-distance physics contained in the form factor f+(s). The B — > Kl + l~ 
matrix elements of four-quark operators, such as (sb)v-A(cc)v-A, are more complicated, 
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but still systematically calculable. Schematically the B — > Kl + l rate is proportional 
to 

ICg + A 49 | 2 + ICxol 2 (35) 

where A 4q represents contributions from four-quark operators, for instance charm loops 
or annihilation effects. In the present discussion we ignore the small contribution from 
CV, which has already been discussed in section 3.2. Typical values are Cg = 4.2, 
Cio = —4.2 and, for the charm-loop amplitude A 4(J m (Ci + 3C2)h(z,s) ~ 0.3. The 
last figure corresponds to an average within < s < 0.25. For large s the charm loop 
develops an imaginary part, but the magnitude is of similar size. Annihilation effects 
are negligible as shown in section 3.3. Thus A 4(J is only about a 10% effect, both as a 
correction to the Cg amplitude and to the total rate. Because this term is numerically 
subleading the impact of any uncertainties in its evaluation will be suppressed. We 
briefly discuss the theory of A 4(J in the regions of low and high q 2 . 

In the low-q 2 region A 4(? can be computed using QCD factorization [28]. This ap- 
proach, which is based on the heavy-quark limit and the large energy of the recoiling 
kaon, should work well for the real part of the amplitude in view of the experience 
from two-body hadronic B decays [32] and B — > K*^y [33]. Power corrections of order 
A/rrib ~ 0.1 in A 4(3 give only percent level corrections for the differential rate. The 
charm loops receive also corrections of order A 2 /ml [34], which have been estimated at 
the level of a few percent for the exclusive decay B — > K*^ [35]. Since the charm loops 
are relatively less important in B — > Kl + l~ by about a factor of five in the rate, the 
impact of the correction is reduced. On the other hand, the effect increases somewhat 
as q 2 approaches the resonance region. In the inclusive case b — > sl + l~ it amounts to a 
few percent [36], averaged over the low-g 2 region. We therefore conclude that the A 2 /m 2 c 
correction is unlikely to affect B — > Kl + l~ in an appreciable way. 

Any quark-level calculation of physical amplitudes involves the concept of quark- 
hadron duality. There are no indications that this assumption, applied to the charm- 
loops for small q 2 up to about 7GeV 2 (s = 0.25), would introduce an error in excess of 
power corrections or perturbative uncertainties [31]. 

Light-quark loops are generally suppressed by small Wilson coefficients (QCD pen- 
guins) or small CKM factors. Violations of local quark-hadron duality could come from 
the presence of light vector resonances at low q 2 . To get an order-of- magnitude estimate 
we consider the branching ratio of the decay chain B~ — > K~p° — > K~e + e~, which is 
measured to be [6] 

B(B~ -> K-p°) x B(p° -> e + e - ) = 

(4.2 ± 0.5) • 10~ 6 x (4.71 ± 0.05) • 10~ 5 = (2.0 ± 0.2) • 10~ 10 (36) 

This contribution arises from the |A 4(? | 2 term in (35) for which quark-hadron duality 
cannot be expected to hold [31]. However, like similar contributions with other light 
vector resonances, it clearly gives a negligible contribution. Resonance effects could be 
more important in the interference of A 4(J with Cg in (35). In this case they are part 
of the hadronic amplitude that is dual to light-quark loops in the partonic calculation. 
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Oscillations in s due to light resonances are bound to be small because of the smallness 
of the light-quark loops. Integration over the low-g 2 region will further reduce such 
violations of local duality to a negligible level. 

In the high-q 2 region the appropriate theoretical framework for the computation of 
A 4q is an operator product expansion exploiting the presence of the large scale q 2 ~ m 2 . 
This concept has been used in [22] in analyzing the endpoint region of b — > sl + l~, 
which is governed by few-body exclusive modes. A systematic treatment, including 
the discussion of subleading corrections, has been given in [37]. Power corrections are 
generally smaller than for low q 2 . Terms of order A/rrib arise at order a s [37] and the 
analogue of the A 2 /m 2 corrections at small q 2 now contribute only at order A 2 jm\ [36]. 
More important are perturbative corrections to the leading-power term, which however 
can be systematically improved. Finally, uncertainties could come from violations of local 
quark-hadron duality. By duality violation we mean deviations of the OPE calculation, 
at fixed q 2 and in principle including all perturbative and power corrections, from the 
real- world hadronic result. Such violations are related in particular to oscillations of A 4q 
in s due to higher charmonium resonances. These oscillations are absent in the smooth 
OPE result. To first order in A 4q only its real part contributes to (35). Implementing the 
higher charmonium resonances, -0(3770), -0(4040), -0(4160), -0(4415), and the hadronic 
cc continuum in the approximation of Kriiger and Sehgal [38], we estimate the relative 
amplitude of oscillations in Cg + ReA^ to be of order 10 to 20%. We expect these 
local variations to be averaged out when the spectrum is integrated over s [22] such 
that the residual uncertainty will be somewhat reduced. The s-integration, which is 
also phenomenologically motivated, effectively produces a smearing that leads to a more 
'globally' defined quantity where duality is better fulfilled. As discussed in [31], global 
duality in this sense cannot be expected to hold for the second order term |A 4(? | 2 in 
(35). On the other hand, this contribution is numerically very small, at the level of few 
percent, and duality violations will only have a minor effect. To illustrate this point 
we consider the decay chain B~ — > i^"-0(377O) — > K~~e + e~, which can be viewed as a 
resonance contribution to |A 4(J | 2 . In the case of the narrow charmonium states a similar 
contribution leads to the very large resonance background mentioned above. Here we 
have [6] 

B(B~ -> #-0(3770)) x 5(^(3770) -> e + e~) = 

(4.9 ± 1.3) • 10~ 4 x (9.7 ± 0.7) • 10~ 6 = (4.8 ± 1.3) • 10~ 9 (37) 

This indicates that resonance contributions are rather small, in agreement with our pre- 
vious remarks. In conclusion, we have argued that duality violations from the resonance 
region at high q 2 are at a moderate level and should not spoil a precision of theoretical 
predictions for (partially) integrated branching ratios of B — > Kl + l~ at the level of sev- 
eral percent. A more detailed investigation of this issue would be of interest and will be 
given elsewhere. 

In the present analysis we ignore higher order electroweak and QED radiative cor- 
rections. The latter could modify the decay modes and their ratios presumably at the 
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Table 3: Input parameters. 



level of several percent. The leading effects could be taken into account if it should be 
required by the experimental precision. 

3.5 D — > K~uu: Background from B — > t~v t — > K.~v T v T 

The decay B~ — > t~v t followed by r~ — > K~v T produces a background for the short- 
distance reaction B~ — > K~vv, which has been discussed recently in [39]. The branching 
fractions of B~ — > t~u t and t~ — > K~v T are given by 

R / R - — \ G 2 F m B m 2 T fl 2 / m 2 \ 2 

5(5 ^rz/ r )=r B Kb 1 2" 

8tt V m Bj 

B(r-^A--^) = r T ^^i|K,| 2 (l-^y 2 = 7.46.10- 3 (39) 

16-7T \ mj J 

The numerical values are based on the input parameters in Table 3. The background 
from the decay chain B~ — > t~v t — > K~v T v T gives a contribution to the dilepton-mass 
spectrum, which can be written as 

~+ ^rgrW R / R - _ - F/ s Ri - _ \ 2t((l-t)(t-r K )-ts) 

(40) 

where we used (4) and t = m 2 ./m 2 B . The spectrum in (40) extends from s = to 
s = (1 — £)(1 — rx/t) = 0.818. The maximum s happens to coincide almost exactly 
with the endpoint of the spectrum in the short- distance decay B~ — > K~uu, s m = 0.821 
[39]. Integrated over the full range in s, the phase-space factor in (40) gives 1. This 
reproduces the result for B(B~ — > K~v T V T \k gr in the narrow- width approximation for 
the intermediate r lepton. 

If the decay sequence B~ — > t~v t — > K~v T v r cannot be distinguished experimentally 
from the short-distance decay B~ — > K^uu, this background should be subtracted from 
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the measured rate of B — > K + "invisible" to obtain the true short-distance branching 
fraction. In the standard model we find from (38) and (39) 

B(B~ -> K-v T v T ) hkgr = (0.65 ± 0.16) • 10~ 6 (41) 

assuming an uncertainty of 25% due to and \V u b\. The central value amounts to 
about 15% of the short-distance branching fraction (44). A subtraction of (41) from the 
measured branching ratio would then lead to an uncertainty of about 0.16/4.4 = 4% 
on B(B~ — > K~vv). This error might be further reduced in the future with improved 
determinations of /# and \V u b\. 

Probably the best method to control the background is to use the experimental 
measurement of B(B~ — > t~u t ). In this way, any new physics component in the latter 
decay will automatically be removed from the measurement of B~ — > K~vv. This will 
simplify the new physics interpretation of the measured B{B~ — > K~vv). The present 
experimental value for B(B~ — > t~u t ) is [7,40,41] 

B(B- -> T-u T ) exp = (1.43 ± 0.37) • 10~ 4 (42) 

Together with (39) this gives 

B(B~ -> K-v T V T ) bkgr = (1.1 ± 0.28) • 10~ 6 (43) 

The 26% uncertainty in (42) thus implies an error of about 6% in B(B~ — > K~vis), 
assuming the central value of (44). However, by the time when £?(£?~ — > K~vv) will 
be measured at a Super Flavour Factory, B{B~ — > t~v t ) will be simultaneously known 
with high precision. According to [3,4] the expected accuracy is about 4% or better. 
Assuming the central values of (42) and (44) as before, the background subtraction will 
then lead to an error of only about 1% in B(B~ — > K~vv). 

We conclude that the background from B~ — > t~u t — > K~u T u T pointed out in [39] 
has to be taken into account for a precise measurement of the short- distance branching 
fraction B(B~ — > K~vv). It needs to be subtracted from the experimental signal, but 
this should ultimately be possible with essentially negligible uncertainty. The background 
discussed here is absent in the case of the neutral mode B° — * K Q vv. 

4 Precision observables 

4.1 Theory expectations for branching fractions 

The input parameters we will use in the present analysis are collected in Table 3. 

To begin our discussion of numerial results we consider first the integrated branching 
ratios of B~ — > K~vv and B~ — > K~l + l~ . For the neutrino mode we find 

b(b- - k-uu) ■ io 6 = 4.4 t\i (um 18:? (<*>) 18:? (&i) (44) 



14 



We have displayed the sensitivity to the form factor parameters, which are by far the 
dominant sources of uncertainty. The form factor normalization /+(0) has the largest 
impact, while the shape parameters are relatively less important. 

The fully integrated, non-resonant B — > Kl + l~~ branching fraction can be evaluated 
in a similar way. This quantity corresponds essentially to the experimental result in 
(3), which has been obtained by cutting out the large background from the two narrow 
charmonium resonances and by extrapolating the measurements to the entire q 2 range 
to recover the total non-resonant rate. The precise correspondence between theoretical 
and experimental results will depend on the details of the cuts and the extrapolation 
procedure. We will treat the resonance region more carefully later when we study pre- 
cision observables. For our present discussion we simply identify the integral over the 
non-resonant spectrum in (8) with the measurement in (3). This appears justified as 
the error from this identification is expected to be below the experimental uncertainty. 
Adopting these considerations we compute 

b(b- -> K-i+r) ■ io 6 = o.58 + _° :ll (/ + (o)) ±8:Sg M ±8:S (h) -om M (45) 

In addition to the still dominant dependence on the form factor we have in this case a 
non-negligible perturbative uncertainty, which we estimate in the standard way through 
a variation of the scale \i between m b /2 and 2rrib around the reference value of /i — m b . 
The scale dependence is at a rather moderate level of ±6% with NLO accuracy, much 
smaller than the error from the hadronic parameters. Within sizeable, mainly theoretical 
uncertainties, the prediction (45) is in agreement with the measurement in (3). 

Whereas the individual branching fractions (44) and (45) suffer from large hadronic 
uncertainties, we expect their ratio to be under much better theoretical control. It 
is obvious that the form factor normalization /+(0) cancels in this ratio. Moreover, 
as illustrated in Fig. 2, the shape of the q 2 spectrum is almost identical for the two 
modes. This is because the additional g 2 -dependence from charm loops in B — > Kl + l~ , 
compared to B — > Kuu, is numerically only a small effect outside the region of the 
narrow charmonium states. As a consequence, also the dependence on the form factor 
shape will be greatly reduced in the ratio 

R ~ B[B- -> K-M-) (46) 

Numerically we find 

i2 = 7.59lg;8}(ao)18:S(6i)SS(|i) (47) 

This prediction is independent of form factor uncertainties for all practical purposes. 
It is limited essentially by the perturbative uncertainty at NLO of ±6%. Using the 
experimental result in (3), the theory prediction (47), and assuming the validity of the 
standard model, we obtain 

B(B~ -> K~vv) = R ■ B(B~ -> K~l + r) exp = (3.64 ± 0.47) • 10~ 6 (48) 
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Figure 2: The dilepton invariant-mass spectra for B — > Kvv (solid) and B — » Kl + l~ 
(dashed). For easier comparison of the shapes the plotted differential branching fractions, 
dB/ds versus s = q 2 /m 2 B , were normalized by their integral B. The reference values have 
been used for all parameters. 



With an accuracy of ±13%, limited at present by the experimental error, this result is 
currently the most precise estimate of B(B~ — > K~vv). 

Since isospin breaking effects in the decay rates are very small, the branching ratios 
for the decays B° — > K°vv and B° —>■ K°l + l~ are given by the branching ratios for the 
corresponding B~ modes multiplied by a factor of t(B°)/t(B~) = 0.93. 



4.2 Precision observables: ratios of branching fractions 

In order to obtain theoretically clean observables, the region of the two narrow char- 
monium resonances ip{lS) and ip{2S) has to be removed from the q 2 spectrum of 
B — > Kl + l~ . This leaves two regions of interest, the low-s region below the resonances, 
and the high-s region above. For the present analysis we define these ranges as 

lows: < s < 0.25 . , 

- - (49) 
high s : 0.6 < s < s m 

The resonance region 0.25 < s < 0.6 corresponds to the q 2 range 7GeV 2 < q 2 < 
16.7 GeV 2 . For our standard parameter set the total rate for B — * Kvv or B — * Kl + l 
(non-resonant) is divided among the three regions, low-s, narrow-resonance, high-s, as 
35 : 48 : 17. 

We first concentrate on the low-s region, where B~ — > K~l + l~ can be reliably calcu- 
lated. To ensure an optimal cancellation of the form factor dependence, one may restrict 
also the neutrino mode to the same range in s and define 

C 25 ds dB(B~ -> K-vv)/ds , 

i? 25 = 50 

/° ds dB{B- -> K-in-)/ds 
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This ratio is determined by theory to very high precision. Displaying the sensitivity to 
the shape parameters and the renormalization scale one finds 

R2 5 = 7M-° ™(a )-° ™(b 1 )- il(ri (51) 

The form factor dependence is seen to cancel almost perfectly in i? 2 5- The shape para- 
meters affect this quantity at a level of only 0.5 per mille. One is therefore left with the 
perturbative uncertainty, estimated here at about ±5% at NLO. 

The independence of any form factor uncertainties in R25 comes at the price of using 
only 35% of the full B~ — > K~vv rate. We therefore consider a different ratio, which is 
defined by 

J° m dsdB(B- -> K-vv)jds 
/„ ds dB(B- -> K-l+l~)/ds + J g ds dB(B~ -> K~l+l-)/ds 

In this ratio the fully integrated rate of B~ — > K~vv is divided by the integrated rate 
of B~ — > K~l + l~ with only the narrow-resonance region removed. This ensures use of 
the maximal statistics in both channels. Due to the missing region in B~ — > K~l + l~ 
the dependence on the form factor shape will no longer be eliminated completely, but 
we still expect a reduced dependence. Numerically we obtain, using the same input as 
before, 

i? 256 = 14.60 (oo) ±8lS§ (61) ;o!62 M (53) 

This estimate shows that the uncertainty from a and b\ is indeed very small, at a level 
of about ±3%. With better empirical information on the shape of the spectrum this 
could be further improved. 

We conclude that ratios such as those in (50) and (52), or similar quantities with 
modified cuts, are theoretically very well under control. They are therefore ideally suited 
for testing the standard model with high precision. 



4.3 Precision observables: B — > Kl + l with lattice input 

Until now our strategy has been to achieve accurate predictions by eliminating the 
form factor dependence altogether. A variant of our analysis consists in taking a single 
hadronic parameter, the form factor at one particular value of q 2 , / + (sq), s o — 9o/ m B> as 
additional theory input. The shape of the form factor can be fitted to the experimental 
spectrum as discussed in sec. 3.2. At the expense of one extra hadronic parameter it is 
then possible to probe short distance physics based on B — > Kl + l~ alone. The necessary 
input /+(so) could come from lattice QCD calculations. 

This approach is analogous to the method pursued in [42] to determine \V u b\ from 
B — > nlis. In this case lattice results on the B — > 71 form factor at a typical value of 
q 2 = 16 GeV 2 were considered as theory input. Experimental data on the spectrum 
and decay rate of B — > ixlv can then be used to extract \V u b\. In order to describe 
the form factor shape, [42] employed dispersive bounds and a related class of general 
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parametrizations [43,44,45]. These more sophisticated parametrizations may also be 
applied in our case, if more than two shape parameters should be required to fit the 
data with the appropriate precision. For the time being the form factor parametrization 
used here is completely sufficient. As shown in [42], the limiting factor is the value of 
/+(so)- We remark that in the case of B — > Kl + l~ , the narrow-resonance region should 
be removed from the analysis when performing the fit to the form factor shape. 

To illustrate the method we extract the form factor f + (s ) at point s = 16 GeV 2 /m 2 B 
from the measurement in (3). Using the best-fit shape parameters a — 1.6 and 61/60 — 1 
we obtain 

U(s ) = 1.05 ± 0.06 (BR) ± 0.03 (54) 

where the first error is from the measured branching ratio in (3) and the second error 
is from scale dependence. The sensitivity to the shape parameters is comparable to the 
uncertainty from the branching ratio. 

Unlike the parameters a and 61, the value of /+(so) determined in this way, assuming 
the standard model, is sensitive to the normalization of the branching ratio and there- 
fore to new physics effects. These could be detected through a comparison with QCD 
calculations of /+(so)- Our default choice of hadronic parameters leads to /+(so) — 116 
but the uncertainty is larger than 15%. Within the coming five to ten years a precision 
of ±4% might be achieved for the form factor /+(so) i n lattice QCD [26]. 

5 New physics 

The branching fractions of B — > Kvv and B — > Kl + l~ are sensitive to physics beyond 
the standard model. If the new physics would modify both of them by (almost) the same 
factor, this change would not be visible in the ratios R25 or i?256 discussed in sec. 4. In 
that case, the new physics could still be seen by studying B — > Kl + l~ separately with 
the method described in sec. 4.3. An example is a scenario with modified Z-penguin 
contributions [5] interfering constructively with the standard model terms. In this case 
i? 25 is changed only by a small amount. 

In general, however, nonstandard dynamics will have a different impact on B — > Kvv 
and B — > Kl + l~ . The excellent theoretical control over the ratios R25 or i?256 will help 
to reveal even moderate deviations from standard model expectations. 

One example is the scenario with modified Z-penguin contributions [5] mentioned 
before, if these contributions interfere destructively with those of the standard model. 
In that case the ratios R25 or i?256 could be significantly suppressed. The modified 
Z-penguin scenario may be realized, for instance, in supersymmetric models [5,46]. 

Another class of theories that do change the ratios are those where B — > Kl + l~ 
remains standard model like while B — > Kvv receives an enhancement (or a suppression). 
Substantial enhancements of BiB — > Kvv) are still allowed by experiment, in fact much 
more than for B — > Kl + l~ . 

A first example are scenarios with light invisible scalars S contributing to B — > KSS 
[46]. This channel adds to B — > Kvv, which is measured as B — > K + invisible. If the 
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scalars have nonzero mass, B — > KSS could be distinguished from B — > fTz/z/ through 
the missing-mass spectrum. On the other hand, if the mass of S is small, or the resolution 
of the spectrum is not good enough, a discrimination of the channels may be difficult. 
The corresponding increase in B{B — > Kvv) could be cleanly identified through the 
ratios R25 and i?256- 

A second example is given by topcolor assisted technicolor [47]. A typical scenario 
involves new strong dynamics, together with extra Z' bosons, which distinguishes the 
third generation from the remaining two. The resulting flavour-changing neutral cur- 
rents at tree level may then predominantly lead to transitions between third-generation 
fermions such as b — > sv T v T . An enhancement of B(B — > Kvv) would result and might 
in principle saturate the experimental bound (1). An enhancement of 20%, which should 
still be detectable, would probe a Z'-boson mass of typically M z > ~ 3 TeV. A similar 
pattern of enhanced B — > Kvv and SM like B Kl + l~ is also possible in generic Z' 
models [46]. 

A more detailed exploration of new physics in B — > Kvv and B — > Kl + l~ is beyond 
the scope of this article. The examples mentioned above illustrate that the ratios of 
branching fractions considered here exhibit a significant sensitivity to interesting new 
physics scenarios. The subject of new physics in b — > svv transitions has been discussed 
in [48] and most recently in [46]. New physics in B — > Kl + l~ has been studied in [49], 
including the information from angular distributions. 

6 Conclusions 

In this paper we have studied precision tests of the standard model through a combined 
analysis of B — > Kvv and B — > Kl + l~ . The main points can be summarized as follows: 

• After removing the narrow-resonance region the branching fraction of B — > Kl + l~ 
can be reliably computed. The dominant amplitude from semileptonic operators is 
simply a calculable expression times the form factor / + . QCD factorization for low 
q 2 and OPE for high q 2 allow one to treat also matrix elements of 4-quark operators 
in a systematic way. These are dominated by charm loops, which are numerically 
small contributions to begin with. Since the tensor form factor fa can be related to 
/ + in the heavy-quark limit, the entire B — > Kl + l~ amplitude becomes calculable 
in terms of practically a single hadronic quantity, the form factor f+(s). 

• The decay mode B — > Kvv is a particularly clean process. It is completely de- 
termined by short-distance physics at the weak scale up to the same form factor 
f+(s). For the charged mode the background due to B~ — > t~v t — > K~v T v T should 
be subtracted from the experimental signal, but this will be possible without in- 
troducing any appreciable uncertainty. 

• The form factor uncertainty can be eliminated by constructing suitable ratios of 
(partially) integrated rates such as i? 25 in (50) and i? 2 56 in (52). The resulting 
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quantities can be computed with high accuracy. The cancellation of form factors 
is exact and does not require the use of approximate flavour symmetries. 

• The perturbative uncertainty of the ratios is estimated to be ±5% at next-to- 
leading order (NLO). This can be further improved by a NNLO analysis along the 
lines of [28]. Uncertainties from other sources are at the level of several percent. 
Some refinements in controlling them should still be possible. 

• Based on the current measurements of B~ — > K~l + l~ we predict 

B(B~ -> K'vv) = (3.64 ± 0.47) • 10~ 6 , (55) 
at present the most accurate estimate of this quantity. 

• The ratios of B — > Kvv and B — > Kl + l~ rates have an interesting sensitivity 
to new physics. The new physics reach benefits from the high accuracy of the 
standard model predictions. A complementary new physics test is possible based 
on B — > Kl + l~ alone, exploiting experimental information on the q 2 spectrum, if 
/+ at one particular value of q 2 is used as input from lattice QCD. 

• New physics in the Wilson coefficients factorizes from low-energy hadronic physics 
in a simple way. The latter is essentially contained only in /+. Our analysis can 
thus be generalized to specific new physics scenarios in a straightforward manner. 

Our proposal puts B — > Kvv as a new physics probe in the same class as K — > irvv, the 
'golden modes' of kaon physics. B — > Kvv together with B — > Kl + l~ thus hold exciting 
opportunities for B physics at a Super Flavour Factory. 



A Relation between form factors fa and / + 

The three form factors in (13) and (14), fr(s), f+(s), and fo(s) or equivalently 

/_( 8 )e[/ ( 8 )-/ + (,)]^ (56) 

q z 

are related in the heavy-quark limit. If we multiply (13) and (14) by ty, = p^/ms and 
use the equation of motion for the heavy quark, fib = b, we find 

2 Z B h = /+-/- (57) 
frig + m K 

Similarly, multiplying (13) with using fib = b, and comparing the result with q^- (13), 
where the quark equations of motion are used on the left-hand side, one finds 

/+ = -/- (58) 
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Together (57) and (58) imply 

Ms) m B + m K 
f+(s) rn B 

The relations (57) - (59) have been obtained in [50] in the heavy-quark limit. They 
apply immediately to the case where the kaon is soft, since then the heavy-quark mass is 
the only large energy scale in the problem. On the other hand, the derivation makes no 
explicit reference to the kaon energy and it has been argued in [51] that these relations 
should be valid in the entire kinematic domain. This conjecture can be justified within 
the framework of soft-collinear effective theory (SCET) [52,53], using the form factor 
relations in the large recoil limit [16,17], which also lead to (57) - (59). A related 
discussion can be found in [54]. 

We may thus use (59) in the entire range of q 2 between and (m B — m K ) 2 . Since 
this expression for j T j f + is an asymptotic result in the heavy-quark limit m b ^> Aqcd, 
an important issue is the question of subleading terms. These can be power corrections 
in Aqcd/^b and perturbative QCD corrections. The perturbative corrections were 
computed in the heavy-quark and large recoil energy limit in [17] with the result (in the 
NDR scheme with MS subtraction) 



m B f T (s) a s (n)C F 



1 



m B + m K f + (s) ' 4tt 



4E K 2E K 
In — k H In 

m b m B — %E K m B 



a s (fi h )C F 4:7r 2 f B f K f 1 ^ 4> K {u) ^ 



4tt Nf + (s)E K X B J l-u 

where fi = (9(m&), fi^ = 0(\/ Ami,), and Ex = 0{m b ) depends on s = q 2 /m 2 B through 
(5). We note that this relation remains valid when the kaon is soft, with Ek = 0{A). 
In that case expression (60) simplifies. The term with a s (nh) is no longer perturbative 
since the effective scale Hh becomes soft. However, the entire term is power suppressed 
~ A/m b because f B ~ l/^/m^ and /+(1) ~ ^ffnb- The second term in square brackets 
is also power suppressed and only the first correction ~ a s \n{fi/m b ) survives. The a s 
corrections can be consistently taken into account at NNLO even though at present, for 
low q 2 , the second term from hard spectator interactions still introduces an uncertainty 
of about 5-10%. 

Power corrections to the heavy-quark limit are more difficult to compute. An estimate 
can be obtained from light-cone QCD sum rules [23], which indicate that fr/ f+ deviates 
from 1 +mjf/mg by less than ±5% for < s < 0.5. The sum rule calculations include 
a s corrections within their framework. 

From these considerations we conclude that the relation (59) should be correct to 
within ±10%. 



B Weak annihilation in B — > Kl + l 

Weak annihilation contributes to B° — ■> K°l + l~ through QCD penguin operators. These 
induce the transition bd — > sd where the valence quarks bd of the B° meson are anni- 
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hilated and transformed into the constituents of the final state K°. The virtual photon 
producing the lepton pair may be emitted from any of the four quarks in this transition. 
In a similar manner the process bu — > su gives rise to weak annihilation in B~ — > K~l + lr , 
where the transition comes from QCD penguins and from doubly Cabibbo suppressed 
tree operators. 

To be specific we treat the case of B° — > K Q l + l~ first. Here the leading-power 
contribution to weak annihilation comes from the QCD penguin operators 

Q 3 = (dtb^v-^Sjd^v-A + ... 

Q 4 = (db) v „ A (sd)v_ A + ... (61) 

The ellipsis refers to similar terms with d replaced by u, c, s and b, which do not 
contribute at the order we are considering. Colour indices are denoted by i, j. 

The kinematics of the annihilation process is conveniently described in terms of two 
lightlike four- vectors n±. Their components can be chosen, without loss of generality, as 

< = (1,0,0,±1) (62) 

The momenta of the B meson, the kaon and the lepton pair, p, k and q, respectively, 
can then be written as 

m B m 2 R — q 2 q 2 tub 

p=— (n + + n_), k = ~^^ q= 2^- B n+ + — n ~ (63) 

For now we assume that the dilepton mass q 2 counts as order Amj, appropriate for the 
low-g 2 region. This means that momentum q is nearly lightlike and approximately in the 
direction of n_, whereas k is lightlike and has the direction of n + if we neglect the kaon 
mass. 

Consider next the bd — > sd annihilation diagram where the virtual photon is emitted 
from the d in the initial state. We will denote by A^i the contribution of this diagram to 
the B° — > K°l + l~ matrix element of Q 4 . Viewed as a function of the d four-momentum 
I = O(A), the diagram has the form 

F(l)=FW(l + ) + r ± FW(l + ) (64) 

up to terms with a relative power suppression in A/m^. The momentum / is decomposed 
into light-cone coordinates l± = n T -l and l±, n± ■ l± — 0, with respect to the two vectors 
n± in (62). The expression for the light-cone projector of the B meson in momentum 
space, appropriate for an amplitude of the type shown in (64), has been derived in [17]. 
It is given by 



,t .fBm El- 
ba = % 



d 



75 (65) 



The derivative in (65) extracts the contribution in (64). After it has been applied, 
l± has to be set to zero and /+ is identified with uj. 



22 



Calculating the contribution A^i with the projector in (65) we find 
An = e 2 f B fK (ujiv) 

Here uj^v is the lepton current in momentum space and Qdi = —1/3 the down-quark 
charge. The second term in (66) is of the same order in A/m& as the first term. It has 
a pole in q 2 and is inconsistent with electromagnetic gauge invariance. However, it is 
structure independent (it depends only on f B , fx, not on the distribution amplitudes) 
and is cancelled by corresponding contributions from the remaining three diagrams. 
These diagrams, where the photon is emitted from the b quark, the s quark and the 
final-state d quark (d2) give explicitly 

A + A + A d2 = -e 2 f B f K Wv) Qb ~ Qs 2 + Qd2 (67) 

q 2 

as contributions to the matrix element of Q4. Charge conservation implies that Qb — 
Q s + Qd2 = Qdi, which guarantees the cancellation of the 1/q 2 term in (66) by (67). The 
first term in (66) then leads to the result (29) once the matrix element of Q 3 is included, 
which is the same as the one of Q<± up to a factor 1/3 from colour. 

We remark that the 1/q 2 term and half of the first term in (66) come from the 
part in (64). The latter contribution, as well as all 1/q 2 terms are absent at leading 
power for weak annihilation in B — > K*j [28,55], in contrast to the present case. 

The result in (29) develops a logarithmic singularity when q 2 becomes soft of order 
A 2 , corresponding to the formal limit q 2 — ■> 0. Since the singularity is integrable, the 
soft region is power suppressed in comparison to the branching ratio integrated from 
to q 2 ~ Ami,, as has been discussed in [28]. 

The annihilation contribution in (32) from operators Q5 = —2(dibj)s-p(sjdi)s+p+- ■ ■ 
and Q e = —2(db)s-p(sd)s+p + ■ ■ • is obtained in analogy to the case of Q 3 , Q4. Also 
for Q 5 , Qq there are structure-independent terms ~ 1/q 2 , which cancel when all four 
diagrams are added. The remaining result is again due to the diagram where the photon 
is emitted from the d quark in the initial state. 
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